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Energy-transfer rate in a double-quantum- well system due to Coulomb coupling 

R. T. Senger and B. Tanatar 
Department of Physics, Bilkent University, Bilkent, 06533 Ankara, Turkey 

We study the energy-transfer rate for electrons in a double-quantum-well structure, where the 
layers are coupled through screened Coulomb interactions. The energy-transfer rate between the 
layers (similar to the Coulomb drag effect in which the momentum transfer rate is considered) is 
calculated as functions of electron densities, interlayer spacing, the temperature difference of the 
2DEGs, and the electron drift velocity in the drive layer. We employ the full wave vector and 
frequency dependent random-phase approximation at finite temperature to describe the effective 
interlayer Coulomb interaction. We find that the collective modes (plasmons) of the system play a 
dominant role in the energy transfer rates. The contribution of optical phonons to the transfer rates 
through the phonon mediated Coulomb coupling mechanism has also been considered. 

PACS numbers: 72.10.-d,73.50.Dn, 73.25.Dx, 73.20.Mf 



I. INTRODUCTION 

Coupled quantum-well systems are known to exhibit 
rich and interesting physics, where correlation effects are 
significant .El In particular, the Coulomb drag effect is a 
unique way of probing majreJbody correlations through 
a transport measurement pu where one of the layers 
is driven by an external current, and the influences 
on the other (drag) layer are investigated. The inter- 
layer carrier- carrier interactions lead to measurable ef- 
fects, such as transresistivity due to momentum transfer 
between the layers. The observed transresistance cru- 
cially depends on the single-particle and collective ex- 
citations of the coupled system, because these excita- 
tions are the ones which mediate the momentum and 
energy transfer between the layers. Thcr|Q-,has been a 
growing theoreticalETH and experimentaloila activity in 
the past years touching upon various aspects of the drag 
phenomenon. 

In this paper we study the energy transfer between 
two layers of quasi-two-dimensional electron gases un- 
der experimental conditions similar to the transresistiv- 
ity measurements. The importance of the energy trans- 
fer between two Coulomb coupled quantum wells were 
pointed out by Price .Ejo In the hot-electron context, the 
energy transfer occurs when there is a difference of elec- 
tron temperatures in the two layers. In the actual drag 
experimentaj'ta the energy transfer was detected from 
the heating effects. The energy transfer rate in spatially 
separated systems wre theoretically considered-also by 
Jacobini and—Price,Ej Laikhtman and .SolomonjLZl Boiko 
and Sirenko,tL3 and recently by Tanatanlj who considered 
the case of a coupled quantum wire system. 

We calculate the temperature dependence of the en- 
ergy transfer rate in a double quantum well system. It 
is assumed that the wells may be kept at different car- 
rier temperatures which are also different from the lattice 



temperature in general.0 The calculations are based on 
the random-phase approximation (RPA), with full con- 
sideration of wave vector and frequency dependencies 
at finite temperatures. The layers are coupled through 
Coulomb interactions, and in the steady state the re- 
sulting charge polarization produces an electrostatic field 
which compensates the drag force in the drag layer. Us- 
ing the momentum and energy balance equations we in- 
vestigate the static and dynamic screening effects on the 
power transfer between the layers. In the drive layer the 
influence of the externally applied electric field is treated 
in terms of the electron drift velocity. We probe the 
effects of a finite drift velocity to study the nonlinear 
regime of the energy transfer rate. The nonequilibrium 
aspects of frictional drag has-recently been considered by 
Wang and da Cunha Lima,Eil who employed the balance 
equation approach. The amount of transfered energy has 
a direct dependence on the electron layer densities which 
can be different in general. We also calculate the effect 
of density mismatch on the energy transfer rates. As a 
final point to be discussed in this report we consider the 
contribution of optical phonon exchange as an additional 
mechanism of interlayer interaction. 



II. MODEL 

We consider two quantum wells of width w, and center- 
to-center separation of d. The potential barriers are as- 
sumed to be infinite, so that there is no tunneling between 
the layers. The two dimensional electron charge density 
in the first layer, m, is related to the Fermi wave vec- 
tor by ra i = kp/(2w), and Tp is the corresponding Fermi 
temperature of the electron gas in the layer. It is also ap- 
propriate to define the dimensionless electron gas param- 
eter r s — v2/ (kpa* B ), where a* B — eo/(e 2 m*) is the effec- 
tive Bohr radius in the layer material with background 
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dielectric constant eo and electron effective mass m* . For 
GaAs quantum-wells a* B ~ 100 A and experimentally re- 
alized electron densities are of the order of 10 11 cm -2 
which corresponds to r s 1 — 2, and Tp « 40— 100 K. We 
take the charge density in the second layer with reference 
to the drive layer; n 2 — an\, so that the quantities for 

(2°) (2) 

the second layer scale as, k F = ^/akp, r y s = r s /^a, 
T ( F 2) = aTp. 

The transport properties of the double quantum well 
system can be characterized by the electron drift veloc- 
ities Vi and electron gas temperatures T^. One of the 
layers (drive layer) is subject to an electric field in the 
x-direction which drives the electrons with a drift ve- 
locity vi . The other well is kept as an "open circuit" , 
therefore v 2 = 0. The drag experiments are performed 
at low electric fields in the linear regime, so we shall 
take the limit v\ — ► at the end of the calculations. In 
this work our starting point for the calculation of the 
energy transfer rate is the balance-equation approach to 
hot carrier transport which has been successfully applied 
to a variety of situations involving transport phenom- 
ena in semiconductors.^ The resulting momentum and 
energy transfer rate equations have-also-been obtained 
within a variety of other techniques Btffl'Ej With the as- 
sumption that only the lowest subband in each layer is 
occupied, the momentum and energy transfer rate ex- 
pressions due to interlayer Coulomb interactions, derived 
within the balance equation approach to nonlinear elec- 
trical trsBWJsort in low dimensional semiconductors, are 
given byc^riy (h = k B = 1); 



/i 2 (i>i - v 2 ) 



E 



LO 



(III; 



n B 



W 12 {q,Lo) 



LO — UO v2 



T 2 



x Imxi(g, w) Imx2 (-9,^12 - , (1) 



and 



/oo 
— LO | W 12 (q,uj) | 5 
-oo " 



LO 

1 77 



n B 



LO - L0 12 

T 2 



x Imxi(g,w) Imx 2 (-g,Wi2 - lo) (2) 

respectively. In the above, lo\ 2 = q x {v\ — v 2 ), Wi 2 (q,Lo) = 
Vi 2 (q) / w) is the dynamically screened interlayer po- 
tential, Imx(<7, lo) is the imaginary part of the tempera- 
ture dependent 2D susceptibility^ for a single layer, and 
ub(x) — l/(exp(x) — 1) is the Bose distribution function. 
The screening function e(q, lo) for the double- well system 
is written as 

e(q,Lo) = [1 - V 11 (q)xi(q,Lo;T 1 )][l - V 22 (q)x 2 (q,^;T 2 )} 
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define the intra- and interlayer unscreened Coulamb in- 
teractions, and Fij(qw) are the form factorstffl for a 
model of infinite barrier and square wells of width w. 
Note that we have indicated explicitly in Eqs. (1) and 
(2) that the layers or quantum-wells are kept at different 
temperatures. Under drag conditions mentioned above, 
the interlayer resistivity (transresistivity) reads 
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where ]\ 2 is the interlayer momentum transfer rate or 
frictional force. Resistivity expression is further sim- 
plified if we consider layer temperatures to be equal, 
T\ = T 2 , and within the linear regime v% — > 0, yielding 



1 



Pl2 



nin 2 e* 



dvi 



(6) 



v 1= 



The energy transfer rate expression given in Eq. (2) 
resembles the momentum transfer rate expression of 
Eq.(l), except the transferred energy lo appears in the 
integrand, and the difference between the Bose distri- 
bution fjuagtions at different temperatures reduce to the 
familiarfiflti - 1/ sinh 2 (w/2T) when T x approaches T 2 . 
With the sign chosen in Eq. (2), P\ 2 is the amount of 
power transferred to the layer 1 from the layer 2. 



III. RESULTS AND DISCUSSION 

We first evaluate the energy-transfer rate Pi2(0) in 
the linear regime (v\ — 0) for a GaAs system. Even 
in this case, the energy transfer rate is non-zero as long 
as the electron gases are kept at different temperatures. 
The comparison of taking the interlayer potential as ei- 
ther statically or dynamically screened is presented in 
Fig. 1(a). Because the energy transfer rate P\ 2 changes 
sign as T2 is scanned for a fixed T\, we plot IP12I in our 
presentations. In the statically screened interaction we 
use e{q) = [1 - Vnxi(q)}[l - ^22X2(9)] - V? 2 xi(q)x2(q), 
in which the static response functions xi. 2 (g,w = 0) en- 
ter. Both quantum wells are taken of width w = 2a* B , and 
with equal electron densities [n\ = n 2 ). The temperature 
of the first layer is kept at T\ — Tp. For three distinct 
sets of r s and d values, the T2 dependence of IP12I is plot- 
ted. We observe that the inclusion of dynamical screening 
effects yields qualitatively and quantitatively different re- 
sults for the energy transfer rate. At very low T2 both 
types of screening yield close results, but with the in- 
creasing temperature, P\ 2 significantly grows for dynam- 
ical screening, whereas it monotonically decreases for the 
case of static screening. The difference between the two 
approaches is attributed to the contribution of plasmons 
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which is completely missed for the statically screened in- 
teraction. Similar differences between the static and dy- 
namic screening approaches were also found in the mo- 
mentum transfer rate at high temperatures determining 
the transresistivity.13 We notice that the qualitative forms 
of the curves are roughly independent of the separation 
distance and the charge density. When the temperatures 
of both layers are equal the energy transfer vanishes for 
all cases in the linear regime. This is also seen in Fig. 1(b) 
where we take T\ at three different values, and vary T2 
from to Tp. Although experimentally the energy trans- 
fer is expected to be from the drift layer to the drag layer 
(in the linear regime at least), for the sake of generality 
the cases where T 2 > T\ are also included in our plots. 
Note that P\i changes sign when T2 > T% , so that energy 
is transferred from the hot layer to the cold one. When 
we compare the overall behavior of P\i in a double-layer 
system with that in a coupled quantum wire system, we. 
observe that in the latter a pronounced peak structureEa 
around T ~ 0.3 Tp is present. On the other hand, the re- 
sults shown in Figs. 1(a) and 1(b) indicate a rather broad 
enhancement coming from the plasmon excitations. Price 
has estimatedli-3 the energy transfer rate (per electron) 
between coupled quantum wells to be ~ 0.1 — 1 erg/s, for 
typical layer densities of n ~ 10 11 cm~ 2 and d « 100 A. 
In obtaining this estimate, the electron temperature was 
taken as 10 3 K, which is about 10Tp. Our calculations 
are mostly done for layers of finite thickness, larger sep- 
aration distances and at relatively lower temperatures 
around Tp. The results indicate rates of the order of 
10 -2 erg/s, which are much smaller than the estimate 
given by Price.liil 
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FIG. 1. The energy-transfer rate for identical quan- 
tum-wells of width w — 2a B , (a) in the static (dashed) and 
the dynamic (solid) screening approximations as functions of 
temperature T2. The temperature of the drive layer is kept 
constant at Ti = Tf- The couple of curves from top to bot- 
tom are for (r 3 ,d) values of (l,5a B ), (2,5a B ) and (2, 7a B ) 
respectively, (b) for different values of T\. The solid, dashed, 
and short-dashed curves are for Ti/Tf = 1, 0.5, and 0.25, 
respectively. 



& 10- 



10"- 



— ■ — 1 — ' — 1 — 1 — 


t — « — 1 — ' — ; 

(a) ' 


\V 


~-T F ] 














1 

V 


\ 

■ 

1 1 1 1 1 


rV 

\ \ j 

1 




0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 

T, / T F T 2 / T F 

FIG. 2. The energy-transfer rate for identical wells 
in the nonlinear regime when (a) T\ — Tf, (b) 
Ti = 0.5Tf- The solid, dashed, and short-dashed curves are 
for vikp/Ep — 0.5, 1, and 2, respectively, in a w = 2a B and 
d = 5a B double- well system (r s — 2). The thin solid curves 
correspond to the linear regime (vi = 0), included for com- 
parison. 

Next we investigate the energy-transfer rate in the non- 
linear regime (i.e. for non-zero v\). The density response 
function and the Bose distribution function of the drag 
layer are calculated at shifted frequencies to — q x v\. In 
Fig. 2 the energy transfer rate is displayed in this non- 
linear situation when T\ is kept at either Tp or 0.5 Tp. 
We observe that for a finite drift velocity and at very low 
temperatures of the second layer, the amount of trans- 
fered power is a few orders of magnitude larger than 
that for the linear regime. As T2 is increased, IP12I 
starts to decrease rapidly and eventually vanish at a crit- 
ical value, before the temperatures of the layers become 
equal. Larger the drift velocity lower the T2 at which no 
energy is transfered between the layers (eg. for T\ = Tp 
and vi = 2E F /k F , P 12 = at T 2 m 0.6T F ). Beyond that 
point, as a consequence of nonlinearity, for further larger 
values of T2 the drive layer start to absorb power from 
the second layer even when T± > T%. The momentum 
transfer rate between two coupled electron-hole quantum 
wells in the jtupnlinear regime was considered by Cui, Lei, 
and HoringBj They found that the nonlinear effects start 
to become significant at different electric field strengths 
(or equivalently the drift velocity v±) for different tem- 
peratures. In a system of two sets of charged particles 
streaming relative to one another the collective modes 
may undergo instabilities with respect to cha#ge density 
perturbations as studied by Hu and Wilkins.cJ It is con- 
ceivable that under the drag effect conditions, such two- 
stream instabilities mayJae detected for large drift veloc- 
ities. Hu and Flensbergca predicted a significant rise in 
the drag rate just under the instability threshold. Recent 
drag rate calculations of Wang and da Cunha Limacil did 
not explore this phenomenon. In our calculations of the 
energy transfer rate, the shift of vanishing P12 point to 
lower T 2 values is expected to reflect the onset of plasma 
instabilities. In our numerical results shown in Figs. 2(a) 
and 2(b) it appears that two-stream instability limit is 
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not yet reached for the drift velocities chosen. 
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FIG. 3. The energy-transfer rate as functions of the 
ratio of drag/drive layer electron densities nz/ni at dif- 
ferent drag-layer temperatures T%, and for electron density 
n\ — 2 x 10 11 cm" 2 (r 3 = 1.22). The dashed and solid curves 
are for the linear and nonlinear (vi = 2Ep/kp) regimes, re- 
spectively. We take w = 2a * B and d — 5a B . Note that in the 
linear regime P12 = when T\ =Tz, regardless of the density 
ratio. 

The effect of charge density mismatch in the two lay- 
ers on the energy-transfer rate is shown in Fig. 3. We 
keep m = 2 x 10 11 cm" 2 (r s = 1.22) and T x = T F con- 
stant, and vary n%jn\ in the range 0.2 to 2 for different 
Ti values. Both in the linear and nonlinear regimes it is 
observed that Pi 2 is very sensitive to the electron den- 
sity ratio of the layers. For relatively large temperatures 
(T2 ~ Tp) the energy-transfer is considerably larger when 
the densities are equal. In the linear regime (cf, dashed 
curves in Fig. 3), the peak value of P12 shifts to smaller 
values of n 2 jri\ ratio when the temperature of the second 
layer is lowered, for example at T2 = 0.1 Tp the largest 
transfer rates occur for 712 ~ 0.2 ri\. For finite drift ve- 
locities the nonlinear effects show up also in terms of the 
density ratio. The vanishing of power transfer at finite 
temperature differences of the layers is seen at even lower 
T2 if n 2 <n\. In the figure, for v\kpjEp = 2, the direc- 
tion of energy flow is from layer 1 to layer 2 (Pi 2 < 0) 
if T 2 < OATp, and vice versa if T 2 > 0.7T F in the cho- 
sen range of the density ratio. For the cases T2 = 0.5 TV 
or 0.6Tp however, Pi 2 < at larger values of n 2 , and 
changes sign at certain values of the density ratio, as n 2 
is lowered. The feature is due to that a constant temper- 
ature correspods to a larger effective temperature for a 
lower-density electron gas, in units of the Fermi temper- 
ature of the second layer. The momentum transfer rate 
in the drag phenomenon has been known to be very sen- 
sitive to the relatijve|-d|ensities in the spatially separated 
electron systems. crOJEJ We find here that energy transfer 
rate also has a strong dependence on the ratio nijn\. 
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FIG. 4. Comparison of Coulomb and optical phonon me- 
diated Coulomb interaction mechanisms in the energy transfer 
rate. The short- and long-dashed curves are, respectively, for 
Coulomb-only and phonon-only contributions to the energy 
transfer rates. The solid curves represent the simultaneous 
inclusion of both mechanisms. The upper and lower curves of 
each type are for d = 5a B and d = 7a* B respectively, (m = ni , 
r s = 2 and w = 2a* B ) 

Finally, we calculate the effect of optical phonons on 
the energy transfer rate in a double-layer electron system. 
At low temperatures, the acoustic phonons-are known to 
contribute to the momentum transfer rate.E3 Since we are 
more interested in the high temperature regime where 
plasmons play a dominant role in the momentum and 
energy transfer rates, we consider only the dispersion- 
less optical phonons for the phonon mediated interlayer 
electron-electron interactions. More specifically, we re- 
place the interlayer effective electron-electron interaction 
by 



W 12 {q,uj) 



Vi 2 (g) + Di2(g,Lo) 
e T (q,uj) 



(7) 



to account for the Coulomb and LO-phonon mediated 
interactions simultaneously. The phonon term Z2ie_which 
also has a Coulombic nature is calculated to bcEj'EZI 



D l2 {q,u) = V l2 (q) 1 



CO 



(8) 



The above form of the interaction due to the exchange 
of virtual phonons arise when the "bubble" diagrams to 
all orders are considered within the RPA. Here 7" 1 is 
the phenomenological lifetime for phonons. 7 is taken 
finite in our calculations but the results do not have a 
significant dependence on 7 as long as it is small (i.e. for 
GaAs, 7 ~ 0.1 meV, whereas wlo = 36meV). The total 
screening function is modified in a straightforward way 
to include the phonon contribution, 

e T (g,w) = [1 - (Vu + £>ii)xi] [1 ~ (V22 + D 22 ) X2 ] 

- (y 12 + d 12 ) 2 XiX2 . (9) 

Da are the phonon mediated intralayer electron-electron 
interaction terms defined similar to Di 2 . The zeros of 
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the total screening function eriq,^) give the coupled 
plasmon-phonon mode dispersions. The purely Coulomb 
and phonon mediated interaction contributions to the 
amount of transfered energy are not simply additive, due 
to the interference terms coming from the total inter- 
action amplitude, Eq. (7) . It is possible to consider the 
phonon contribution alone by keeping only the -D12 term 
in the numerator of Eq. (7L-,i.e. Wi2(q,u}) — D^/st- 
However one should retained the coupling of Coulomb 
and phonon terms in £t(<Z,w). 

We calculate the energy transfer rate between two 
quantum wells by modifying the expression for Pi 2 as set 
out above for the phonon mediated Coulomb interaction. 
We assume that electrons and phonons are in thermal 
equilibrium at T\ . In Fig. 4 we display the energy-transfer 
rates due to the Coulomb and phonon coupling inter- 
actions alone, and when both mechanisms are present 
together. It is observed that the phonon contribution, 
being almost an order of magnitude smaller, has quali- 
tatively similar form as the direct Coulomb interaction 
contribution. The mentioned interference is observable 
when the temperature difference of the layers is small; the 
combined effect of two mechanisms can lead to slightly 
smaller transfer rates (cf, Fig. 4(a)). 

In summary, we have considered the energy transfer 
rate in a double-quantum-well system in a drag exper- 
iment type setup. The interlayer Coulomb scattering 
mechanism through dynamical screening effects greatly 
enhances the energy transfer rate from one layer to an- 
other. We have found that a large drift velocity corre- 
sponding to large externally applied field greatly modifies 
the energy transfer rate and may lead to power absorp- 
tion from the cooler electron gas to the other. Some of 
our predictions may be tested in hot-electron photolu- 
minescence experiments, in which the power loss of an 
electron gas is measured. Lastly, we have also consid- 
ered the contribution of optical phonons in the phonon 
mediated Coulomb drag which may be important to un- 
derstand the future experiments. It would he-interesting 
to extend the measurements of Noh et alBB to higher 
temperatures to observe the effects of coupled plasmon- 
phonon modes. 
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